Nonlinear Bloch-wave interaction and Bragg scattering in optically- induced lattices 
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We study, both theoretically and experimentally, the Bragg scattering of light in optically-induced 
photonic lattices and reveal the key physical mechanisms which govern nonlinear self-action of narrow 
beams under the combined effects of Bragg scattering and wave diffraction, allowing for selecting 
bands with different effective dispersion. 
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Dispersion and diffraction are the fundamental phe- 
nomena of wave physics, and they are responsible for 
temporal pulse spreading and beam broadening. The 
rate of wave broadening can be controlled in periodic 
structures where the effective geometric dispersion pro- 
vides a key physical mechanism for manipulating waves in 
various physical systems, including the Bragg gratings in 
optical fibers 1], waveguide arrays 0, and Bose-Einstein 
condensates (BECs) in optical lattices Q. 

A strong change of the wave dispersion and diffrac- 
tion occurs in the vicinity of spectrum band gaps where 
waves experience resonant Bragg scattering from a peri- 
odic structure |4j . In the structure shown in Fig.QJ both 
scattering and spreading depend on the angle between 
the beam and grating. For broad beams, the evolution 
of transmitted and scattered waves is defined by the effec- 
tive dispersion parameters at the spectral peaks |5( , and 
can be studied in the effective-mass approximation Q. 
The dynamics of narrow beams and wavepackets becomes 
more complicated, and it is strongly affected by nonlin- 
earity, as recently discussed in Ref. 3]. However, the key 
physical mechanisms which govern nonlinear self-action 
of narrow beams under combined effects of the Bragg scat- 
tering and wave diffraction remain largely unexplored. 

In this Letter we study, both theoretically and experi- 
mentally, light propagation in optically-induced photonic 
lattices near the Bragg scattering angle. For the first time 
to our knowledge, we reveal generic relations between the 
Bloch-wave spectrum of the periodically modulated me- 
dia and the specific structure of the diffracted waves de- 
scribing analytically the key patterns of the beam scatter- 
ing and self-action observed in our experiments. We also 
study the nonlinear Bloch-wave interaction and observe 
self-focusing in the spectrum gaps tracing the signatures 
of the Floquet-Bloch solitons. 

Model. Essential features of light scattering in periodic 
photonic structures and the BEC dynamics in optical lat- 
tices can be described by similar mathematical models. 
For definiteness, we consider the optical beam propaga- 
tion in a one-dimensional periodic structure described by 
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FIG. 1: Schematic of a laser beam propagating in a periodic 
photonic structure. The input Gaussian beam is split into 
several Bloch modes due to the Bragg scattering; top plot 
shows the output profile superimposed on top of the grating. 
Beam evolution is defined by the inclination angle a. 



the dimensionless equation, 
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v(x)E + f(x,\E\ 2 )E = 0. (1) 



where E{x, z) is the normalized electric field envelope, 
x and z are the transverse and propagation coordinates, 
normalized to the characteristic values x s and z 3 , respec- 
tively, D = zo\/(4:7rnoxl) is the beam diffraction coeffi- 
cient, no is the average medium refractive index, and A 
is the vacuum wavelength. Function v(x) characterizes 
the periodic modulation of the optical refractive index, 
v(x + d) = v(x), where d is the spatial lattice period. 
The function T{x, |£^| 2 ) defines the self-induced changes 
of the refractive index. 

The theoretical analysis presented below is equally ap- 
plicable to the BEC dynamics including the recent ex- 
perimental results for the BEC dispersion management 
in optical lattices p|. In this case, E(x,z) stands for 
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FIG. 2: (a) Dispersion curves for the Bloch waves in the 
first and second transmission bands. Shading marks the 
band gaps. (b,c) Propagation angle and diffraction coeffi- 
cients of the Bloch waves, (d) Excitation coefficients of the 
Bloch waves; shading shows the normalized Fourier spectrum 
(\F(k)\ 2 ) of a beam incident at the Bragg angle. All parame- 
ters correspond to the experimental conditions. 



the mean-field condensate wave function, x is the spa- 
tial direction along a cigar-shaped trap, z is the normal- 
ized time, and D is a constant parameter. The periodic 
function v(x) defines the effective potential of an optical 
lattice, and T{x, |£?| 2 ) — 7|^| 2 characterizes the inter- 
atomic scattering. 

General analytical approach. First, we present our gen- 
eral theory for describing scattering and diffraction of 
optical beams in periodic photonic structures [see Fig. E] 
based on the Floquet-Bloch formalism pj. The Bloch 
waves are special solutions of the linearized equation Q) 
of the form, E Kin (x,z) = ip Ki n(x) exp(mx + if3 K , n z), 
where /3 KiU and n are the Bloch- wave propagation con- 
stant and wavenumber, respectively, and the index n = 
1,2, .. . marks the transmission band. The Bloch wave 
^K,n(%) has the periodicity of the photonic structure, 
ipK,n(% + d) = i/j KiTl (x), and this constraint defines a spe- 
cific dispersion relation /3« jn , as sketched in Fig.J^a). In 
a one-dimensional structure, the transmission bands are 
separated by gaps 4] , where the wave propagation is for- 
bidden (k is real). The band-gaps are associated with 
Bragg reflection of light at certain incident angles. 

The Bloch waves form a complete orthogonal set of 
eigenfunctions, and the linear evolution of an input beam 
in a periodic lattice is fully defined by its Bloch-wave 
spectrum B(k,,ti). Indeed, the electric field envelope E 
can be found at a fixed z as follows 

E(x,z) = J2 B(^n)i/j^ n (x)e iKX+i ^ z dK. (2) 

„ J-ir/d 



The Bloch-wave spectrum can be linked to the Fourier 
spectrum of the input beam, 

) = C n(K + 27rra/cf)FO + 27rra/d), (3) 



where F(k) = (27r) _1 E(x, 0) exp(— ikx)dx, and 
C u {k) are the Bloch-wave excitation coefficients, as 
shown in Fig. Efd). 

In order to capture the essential features of the beam 
diffraction in a periodic lattice, we calculate the beam 
profiles in the far field regime. Assuming that z exceeds 
several diffraction lengths, we use the stationary-phase 
approximation [3 to separate the leading-order contri- 
butions in the integral representation (0 and obtain 



E(x, z) ~ 



B(K,n)>ifj K}n (x) 

(I0$l*) 1/2 



(4) 



where <p = kx + (3 Kin z + (7r/4)sign(/3i 2 n), and = 
& pK,n/dtf • Solution allows us to link the spatial pro- 
file of the diffracted beam and its Bloch-wave spectrum, 
and it provides a nontrivial generalization of the well- 
known fact that in a homogeneous medium the profile of 
a diffracted beam maps directly its Fourier spectrum • 
According to the summation condition in Eq. (@J), $t} n 
defines the propagation angle (or the group velocity) of 

the Bloch waves. The second derivative of the Bloch- 

(2) 

wave dispersion, /%,n, characterizes the rate of the beam 
spatial divergence, and it has a meaning of the effective 
diffraction coefficient. In the regions where < 0, 

the beam diffraction is normal (as in a homogeneous 

(2) 

medium), and it is anomalous for /%,n > 0. As illus- 
trated in Figs. [21 a-c), the group velocity vanishes at the 
band edges, and it reaches a maximum absolute value 
inside the band, whereas the beam diffraction changes 
from normal to anomalous. A fundamental conclusion 
of our analysis is the following: at any given spatial lo- 
cation x, the beam profile is found as a superposition of 
Bloch waves from different spectral bands, and each of the 
bands gives two contributions associated with both nor- 
mal and anomalous diffraction. Nonlinearity produces 
the opposite effects on such waves, resulting in either 
self-focusing or self-defocusing, and it is the competition 
of these effects that will define the beam shaping in pe- 
riodic structures. 

Equation Q becomes formally invalid in the vicinity of 
a zero-diffraction point, which corresponds to a diffrac- 
tion catastrophe. The corresponding contribution to the 
total field can be expressed through the Airy functions 

by taking into account higher-order diffraction terms. 

(2) 

When Pk/u = 0, the general expression can be simplified, 
-27rz- 1 /3 jB ( /t , n )^ n ( a; )| 6 /4 3 ) i |i/3[ r( _i/3)]-i. The 

contribution from the zero-diffraction region may become 
dominant in the far field where B{K J i,n)\zf3^'^ n \ 1 ^ <C 
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B(n2^n)\zf3 { k^ n \ 1 ^ 2 . This effect is responsible for the for- 
mation of a specific diffraction pattern with two side lobes 
in waveguide arrays |p. In the case of dynamically in- 
duced gratings 0, the refractive index contrast is 
smaller, and the spectral regions with near-zero diffrac- 
tion are relatively narrow [see Fig. E(c)] and are not ex- 
pected to make a dominant contribution for experimen- 
tally feasible propagation lengths. 

Experiment vs. theory. In order to study the Bragg 
scattering and Bloch-wave interaction experimentally, we 
create an optically-induced periodic lattice by interfer- 
ing two ordinary-polarized coherent laser beams in a 
biased photorefractive crystal SBN:60, as discussed in 
Refs. tJLLII- We apply an external electric field (up to 
4000 V/cm) to the crystal, creating a periodic modula- 
tion of the optical refractive index through the electro- 
optic effect. In the SBN crystal, orthogonally polar- 
ized waves have significantly different electro-optic coef- 
ficients. We use the ordinary polarized interfering beams 
which propagate linearly and create a stable periodic lat- 
tice, whereas the extraordinary polarized probe beam ex- 
periences strong nonlinear self-action, which can be quan- 
tified through the nonlinear change of the refractive in- 
dex, v(x)+F(x, \E\ 2 ) = -j(I b + I p (x) + |£| 2 )"\ where 
lb is the constant dark irradiance, I p (x) = I g cos 2 (7rx/d) 
is the two-beam interference pattern which induces the 
periodic lattice, and 7 is the nonlinear coefficient which is 
proportional to the applied DC field. In order to match 
the conditions of our experimental observations in the 
theory and numerical calculations, we use the follow- 
ing parameters: A = 0.532/mi, no = 2.4, xo = lum, 
zo = 1mm, d = 15.1, lb = 1, I g = 0.4, 7 = 2.13, the 
crystal length is L = 15mm, and the input beam width 
is 21.16/im. A detailed description of our experimental 
setup can be found in our earlier paper [To| . 

First, we perform numerical analysis of beam scat- 
tering in the linear regime using the exact Bloch-wave 
decomposition [Eq. Q]- When the inclination angle of 
the scattered beam approaches that of the Bragg reso- 
nance, we observe the simultaneous excitation of several 
spectral bands by an input Gaussian beam, as illustrated 
in Figs. Efa,b). These results are in perfect agreement 
with our general conclusion that a spatial profile of the 
diffracted beam maps various Bloch-wave spectral com- 
ponents. Indeed, three spectral regions can be excited 
simultaneously according to Fig. Efd), where we super- 
impose the Fourier spectrum of the input beam on top 
of the Bloch-wave excitation coefficients. 

In experiment, first we study the linear scattering when 
the beam power is low. Experimental results are summa- 
rized in Fig. ^ where we show the dependence of the out- 
put beam profile on the tilt angle of the periodic lattice. 
Close to normal incidence, the first band is dominant, 
whereas at larger angles the excitation of the second band 
grows. In the vicinity of the Bragg angle [marked as a 
dashed line in Fig.^Jleft)], the transmitted beam consists 
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FIG. 3: Numerically calculated output beam profiles in the 
linear regime for an input Gaussian beam incident at the an- 
gles (a) 10% below and (b) 10% above the Bragg resonance. 
Intensity levels are normalized to the input peak intensity Jo- 
Shown is decomposition of the total intensity profile (top) be- 
tween spectral regions with normal and anomalous diffraction 
of Bloch waves in the first and second bands (as indicated by 
labels). Shading marks the positions of the grating minima. 
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FIG. 4: Experimentally observed beam scattering by an 
optically-induced lattice for low laser powers. Left: output 
intensity vs. the tilt angle; right: beam profiles at different 
angles, shading marks the positions of the grating minima. 



of two parts which correspond mainly to the excitation 
of the first and second bands, as predicted theoretically 
(cf. Fig. |HJ). For angles larger than the Bragg angle, 
the transmitted beam is composed mostly of the Bloch 
waves from the second spectral band. We note that the 
result of scattering very weakly depends on the initial 
beam position. This happens because the spectral width 
of the input beam is less than the width of the first Bril- 
louin zone (Afc < 2ir/d), and according to Eq. Q the 
beam shift introduces a uniform phase modulation of the 
Bloch-wave spectrum but does not change its shape. 
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(a) Bragg angle - 10% 



(b) Bragg angle + 10% 
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FIG. 5: Numerically calculated output beam profiles for dif- 
ferent values of Iq. Beam parameters correspond to Fig. 



(b) Bragg angle + 10% 
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FIG. 6: Experimentally observed output beam profiles for 
different intensities, (a) below and (b) above the Bragg angle. 



As the beam intensity grows, nonlinearity leads to 
strong changes in the output beam profiles, see results 
of numerical simulations of the full model Eq. in 
Figs. |5Ja,b). In the nonlinear propagation regime, nor- 
mally diffracting Bloch-wave components tend to self- 
focus whereas anomalously diffracting components expe- 
rience self-defocusing. This fundamental property of the 
beam scattering suggests a possibility of performing non- 
linear Bloch-wave spectroscopy of the optically-induced 
photonic lattices. In particular, we find that the beam 
self-action effect depends strongly on the incident an- 



gle. When the angle is below the Bragg resonance [see 
Figs. EJa)], the normally diffracting mode is dominant 
in the first band, see Fig. CJa). Therefore, as the beam 
power grows, the mode experiences self-focusing. At the 
same time, we observe a strong nonlinear ity-induced en- 
ergy transfer into the right part of the beam, which corre- 
sponds to the excitation of the second band. In contrast, 
for the scattering angles above the Bragg resonance [see 
Figs. Efb)], the excitation of the second band is domi- 
nant, and in the first band anomalously-diffracting mode 
is stronger, see Fig.OJb). This leads to defocusing on the 
left part of the transmitted beam, whereas the right part 
self-focuses. At higher powers, the beam starts erasing 
the optically-induced grating, and scattering is signifi- 
cantly reduced for both incident angles. In experiment, 
we increase the intensity of the probe beam, and observe 
the transformations of the beam profiles at the crystal 
output [see Fig. |§J in a remarkable agreement with the 
theoretical predictions [cf. Fig. Ej. Numerical calcula- 
tions indicate that either one or two Floquet-Bloch soli- 
tons 11] can form under appropriate conditions, however 
the soliton formation distance exceeds the crystal length 
due to strong Bloch-wave interactions. 

In conclusion, we have studied the Bragg scatter- 
ing and nonlinear Bloch-wave inter-band interaction in 
optically-induced photonic lattices. We have established 
generic relations between the Bloch-wave spectrum and 
the structure of diffracted waves describing theoretically 
the specific patterns of the beam scattering and self- 
action observed in experiment. 
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by the Australian Research Council. 
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